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Abstract: Light scattering detection in the near field, a rapidly expanding 
family of scattering techniques, has recently proved to be an appropriate 
procedure for performing dynamic measurements. Here we report an inno- 
vative algorithm, based on the evaluation of the Exposure Time Dependent 
Spectrum (ETDS), which makes it possible to measure the fast dynamics 
of a colloidal suspension with the aid of a simple near field scattering 
apparatus and a CCD camera. Our algorithm consists in acquiring static 
spectra in the near field at different exposure times, so that the measured 
decay times are limited only by the exposure time of the camera and not by 
its frame rate. The experimental set-up is based on a modified microscope, 
where the light scattered in the near field is collected by a commercial 
objective, but (unlike in standard microscopes) the light source is a He-Ne 
laser which increases the instrument sensitivity. The apparatus and the 
algorithm have been validated by considering model systems of standard 
spherical nano-particle. 
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1. Introduction 

During the last few decades, light scattering techniques have demonstrated their ability to pro- 
vide detailed information about properties of complex fluids. This capability is mainly due to 
the fact that these methods allow direct measurement of ensemble averages of the spatial and 
temporal fluctuations of fluid properties, which are precisely the quantities that can be usually 
calculated by theoretical methods of statistical mechanics, thus making available a simple way 
to compare experimental results with theory. 

More recently, low angle light scattering methods have also became very popular thanks to 
the advent of pixilated sensors like CCD or CMOS cameras and their continuous technical 
development. The utilization of such sensors was first proposed by Wong and Wiltzius in 1993 
JT] and gave rise to a large number of similar approaches (0 [3J. In these first attempts to 
utilize a multi-element sensor to gain statistical accuracy (typical sensors have about lk x lk 
pixels, thus providing 10 6 independent measurements at one shot), the sensor was positioned 
in the Fourier plane of a lens. In this way it was possible to collect the scattered light in the 
far field, after removing the probe beam by means of a small mirror in the center of the focal 
plane. Accordingly, we will refer to this family of procedures as Scattering In the Far Field 
(SIFF) techniques. In order to overcome some of the difficulties experimented with SIFF, an 
alternative approach was recently proposed, consisting in collecting the light scattered by the 
sample in the near field, i.e. very close to the sample. We will refer to this second family of 
techniques as Scattering In the Near Field (SINF). The SINF methods show some advantages 
with respect to SIFF methods, including the extension of the scattering angle range to arbitrarily 
small values and the avoidance of stray light problems. Examples of SINF techniques are the 
shadowgraph (4][5]|6), the schlieren, the speckle schlieren Q, and the Heterodyne Near Field 
Scattering (HNFS) 0. In the present paper, we will consider heterodyne SINF techniques 
only, in which both the scattered beam and the more intense transmitted beam are conveyed 
to the sensor. In this case, Fourier analysis of the heterodyne SINF images allows to recover 



the intensity of the scattered beams from the interference fringes they generate on the nearly 
uniform transmitted beam. 

SINF procedures have also been extended in the time domain in order to perform Dynamic 
Light Scattering (DLS) measurements (9l LT21 HH LT21 ^ offering detailed data of the time correla- 
tion functions of the non-equilibrium fluctuations or of the diffusional random walk of colloidal 
particles. This approach has also been applied to optical microscopy fl3l . suggesting that, even 
with incoherent illumination, one can achieve information about sample dynamics. 

When using SINF techniques the dynamic data are obtained by acquiring sequences of im- 
ages at a given frame rate. After the acquisition, the images are Fourier transformed and pro- 
cessed. In the literature two different algorithms have been used to analyze the sequences, 
namely: (1) calculation of the time-correlation function (or equivalently the time power spec- 
trum) and (2) calculation of the time structure function. In the present paper we introduce and 
describe a third method, which allows measuring fast dynamics without the need for a fast 
acquisition. 

Approach (1) is fairly similar to the standard procedure used in DLS experiments, where 
the time-correlation function is obtained by hardware calculation using the data collected by 
a photo-multiplier tube aligned at a certain angle with the probe beam |[T4| . For example, this 
approach has been used on shadowgraph images to investigate the dynamics of spiral defect 
chaos in high pressure CO2 lfl5l and also to analyze the propagating modes that appear in a 
fluid subjected to a temperature gradient lfl6l . 

Method (2) originated from the pioneering proposal by Schulz-DuBois & Rehberg ifTTl to 
develop hardware (a so called "structurator") capable of computing the structure function in 
place of the standard "correlator". This technique permits to eliminate the annoying slowly 
moving drifts that would plague scattering experiments. This was first reported by some of us 
J9][l0|, and was later applied to a broad variety of different experimental methods like X-ray 
scattering lfT8l . optical microscopy liT3l and HNFS lTT2l . 

The innovative approach proposed in the present paper is based on the study of the de- 
pendence of the (static) image power spectrum on the camera exposure time, hence the name 
Exposure Time Dependent Spectrum (ETDS). The physical idea behind this technique is that, 
when increasing the exposure time, the fast dynamics of the images is progressively faded out 
and what is left are the long time fluctuations only. In practice, increasing the exposure time 
is equivalent to the application of a sort of low pass filter, averaging the fastest fluctuations 
and thus depleting the power spectrum. For very long exposure times the fluctuations are com- 
pletely smeared and we detect the background electronic noise of the CCD camera. Normally, 
state-of-the-art methods can detect fast dynamics up to the camera frame rate; by contrast, our 
ETDS method works by acquiring uncorrelated images at a very low frame rate, and the fastest 
dynamics it can detect is limited only by the shortest exposure time the camera can reach. This 
methodology intrinsically boosts the performance of any camera, since the shortest gating time 
is always much shorter than the fastest frame-rate delay. We suggest that the method can yield 
very impressive results, down to the sub-picosecond scale, by using pulsed lasers or fast-gated 
cameras or intensifiers. The method has been validated by measuring the time constants of three 
colloidal particles of different sizes using a HNFS set-up implemented on a microscope. 

The rest of the paper is organized as follows: in section 2 the experimental apparatus used to 
collect scattering images will be described; in section 3 the ETDS will be derived; in section 
4 the experimental data for a wide range of wave vectors will be presented; finally, a brief 
summary of the work will be given. 
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Fig. 1 . Scheme of the optical set-up. A plane parallel beam of a He-Ne laser is attenuated by 
a neutral filter and bend upwards by a mirror. The beam is expanded by means of a negative 
focal length lens, making it diverge slightly before it enters the sample cell. Scattered light 
is acquired in the near field together with transmitted light, through a 40 x microscope 
objective, which conjugates a plane close to the sample onto the CCD sensor. 

2. Experimental set up 

In Fig.[T]a sketch of the optical set-up is presented. It is based on a modified microscope where 
the light source is a He-Ne laser (10 mW, NEC) with an output of 0.7 mm diameter TEM00. The 
beam is used directly, without spatial filtering, because the quality of the unmodified beam has 
been proved to be enough for our purpose. The laser beam passes through a neutral filter wheel, 
with transmission range 0.3 - 0.0003. After the filter, the beam is reflected upward by a mirror, 
and goes through a -5 cm focal-length, negative double-concave lens, in order to increase its 
diameter. We have checked that the slight divergence of the light impinging on the sample can 
be neglected. The sample is placed 15 cm after the lens, where the beam diameter is about 
2.5 mm wide. 

The sample is placed in a glass cell with optical path 1 mm, made of two microscope cover 
slips, spaced by small glass strips cut from a microscope slide, and glued with silicone rubber. 
Over the top cover slip, we placed an iris with a 1 mm-diameter aperture, to reduce unwanted re- 
flections. The optical system consists of a plan-achromatic 40x objective (Optika Microscopes, 
FLUOR), with 0.65 numerical aperture, 160 mm focal distance and 0.17 mm working distance; 
its focus lies on plane of the iris, about one mm outside the cell. Images are acquired with a 
CCD camera (Andor Luca), whose sensor is 658 x 496 pixels. The camera's maximum frame 
rate is 30 frames per second, with a minimum exposure time of 0.5 ms. The sensor is placed at 
160 mm from the microscope objective, so that it collects images directly with a magnification 
of about 40x. The CCD sensor images an area of 200 x 150/xm 2 , while the diameter of the 
illuminated area of the sample is 2.5 mm, thus providing the conditions for near field detection 
of the light [12]. Each experiment consists of a set of measurements collected at different expo- 



sures times. For each measurement a different neutral filter is selected, and the exposure time 
is adjusted so that the acquired images have an average intensity corresponding to half of the 
dynamic range. One hundred images are then acquired, at a frame rate of 1 Hz, slowly enough 
to ensure that all the images are uncorrected. 

3. Image processing 

The goal of static scattering measurements is to obtain information about /(Q), the scattered 
intensity at transferred wave vector Q, where Q = k s — k,, k s is the wave vector of the scattered 
beam and k, is the wave vector of the incident beam. In conventional static SIFF experiments, 
the static scattered intensity / (Q) is obtained by direct measurement, collecting the scattered 
light by using a photo-multiplier tube or a CCD in the far field at various angles with the 
impinging beam. By contrast, in a static SINF experiment using a pixilated sensor, we measure 
instantaneously the image signal i (x), which is the intensity distribution at point x on a plane 
close to the sample. After a set of images is acquired, the Image Power Spectrum (IPS) Si (q) 
is evaluated, calculating the average of the square modulus of the Fourier transform. 

5 ; -(q)-(|'(q)| 2 ) = (l^['W]| 2 ) (1) 

where & ( • ) is the 2D Fourier transform, q is the 2D spatial wave vector on the image plane, and 
the mean value (■) is obtained by averaging over the set of images. Usually, before this Fourier 
processing, the images are "cleaned" by subtracting the optical background, which is obtained 
by averaging the images. The IPS S,-(q) represents the intensity of the Fourier modes, each 
one corresponding to a diffraction fringe generated by the interference between the transmitted 
beam and the scattered beam. By calculating the projection of Q on the image plane, it's easy to 
show that the transferred wave vector Q is related to the image wave vector q by the following 
relation: 

Q (q) = V2kJ 1 - \/Hf7 (2) 

where k as lc s as k[ is the light wave vector in the medium. The relation between the IPS S\ (q) 
and the scattered intensity / [Q (q)] = / (q) is linear, and is actually given [9, 1 1 by the sum of 
two terms: 

5 ( -(q)=/(q)-r(q)+B(q) (3) 

where T (q) is the instrument transfer function, i.e. a relation between the wave vector and 
the instrument sensitivity, and B (q) is the electronic background accounting for noise sources 
within the grabbing process. The instrument transfer function varies depending on the exper- 
imental set-up; it can be simply equal to a constant in the case of HNFS (8] QjO or schlieren 
ll7l [T0l thus providing the static scattered intensity with no further complications. Conversely, it 
has been shown [5,6, 1 1 1 that for a shadowgraph the transfer function exhibits, in some q range, 
deep oscillations, thus making the retrieval of scattered intensity difficult and providing no in- 
formation in that vector range. In other cases the instrument transfer function can be completely 
unpredictable, which makes the static data completely unusable. However, the knowledge of the 
transfer function is not needed for the dynamic analysis, in which the attention is focused on 
the time fluctuations of the scattered intensity. 

Indeed, the goal of dynamic light scattering is to study the time correlation function of the 
field Ce (q,f )■ In the case of dynamic SIFF in heterodyne configuration, this quantity is usually 
directly measured with a correlator; or in the case of homodyne configuration, through Siegert's 
relation. Our ETDS-based algorithm for SINF analysis relies on the following argument. As 
pointed out by Oh et al. ll20l . when the camera exposure time Af is not negligible, the IPS 



Sj (q) depends also on At. Hence we will call it ETDS and use the notation S,- (q, Af). It is easy 
to find a relation between the field correlation function Ce (q,?) and Si (q, Af) as follows. An 
image i(x,At), obtained with an exposure time Af, is the results of a time averaging over the 
instantaneous intensity map i (x,f ): i (x, Af) = ^ J/° i (x,?')df', where to is the starting time of 

the image exposition. The same applies for the Fourier transform: i (x, Af ) = ff° +At i(x,t') dt'. 
The ETDS is then formally obtained as follows: 

/ i\ 1 / /"'o+ A ' /"'o+ A ' 

S ; -(q,Af) = (|/(q,Af)| 2 ) = ^^ J dt"i (q/) i (q,t") 

I rtQ+tx rtQ+M 
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where C\ (q,f ) is the time correlation function of the Fourier modes of the images. The double 
integral can be reduced to a simple integral with the following change of variables: 
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Substituting we have: 
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The calculation of the integral over ds" gives: 
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Si(q,At) = -^J ds(At-s)Q(q,s) (7) 



but, in turn, Q (q,f) is related to Ce (q,?) and so: 
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Si (q, Af ) - T (q) ^ ) ^ (Af - s ) C £ (q, s) + B (q) (8) 

This is the general relation between 5, (q, Af ) and the correlation function Ce (q,f ): the ETDS 
calculated at a certain exposure time Af is the average of the field correlation function, weighted 
by a triangular function which vanishes for s — > Af. This implies that spectra obtained with 
different exposure times give different results, and the analysis of this variation as a function 
of Af brings information about the sample dynamics, or the field correlation function Ce (q,f )■ 
Now we consider the specific case we studied in our experiments: colloidal particles performing 
a Brownian motion. In this case the field correlation function is a decreasing exponential: 

C E (q,t)=I(q)exv(- t -) (9) 



The resulting ETDS is: 
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The last equation indicates that S,-(q,Af) is a decaying function with a characteristic time T. 
In general this time constant can be different for different wave vectors Q. For example for 
diffusing Brownian particles the time constant varies as %{Q) = 1/ {DQ 2 }, in which D is the 
translational diffusion coefficient of the particles in the solvent, as given by the Stokes-Einstein 
equation. 

It's worth analyzing the limits of the ETDS for very short and very long exposures times: 

lim 5,(q,Ar)=/(q)r(q)+B(q) (12) 

Af/T<l 

lim Si(q,At)=B(q) (13) 

A//t»1 

For very short exposure times the ETDS approximates the static value of the IPS as given by Eq. 
(0. For long times the ETDS goes asymptotically to the electronic background of the system. 

In summary, by acquiring series of images with many different exposure times and analyzing 
the ETDS for each single wave vector by using the fitting procedure in Eq. dTOb and (fTTT >. one 
obtains a direct and simultaneous measurement of the time constants of the sample for all 
the measured wave vectors. As a byproduct, one gets also two other terms from the fitting 
procedure, namely, the product of the static scattered intensity times the instrument transfer 
function / (q) ■ T (q) and the electronic background of the system. 

4. Experimental results 

In order to test the validity of the our approach we have performed three experiments with 
different commercial samples of polystyrene nano-particles of 80 nm, 150 nm and 400 nm 
diameter (Polyscience Inc.). The particles are well mono-dispersed and the declared size cor- 
responds to the measured size within 5% (81±4 nm, 149±5 nm, 402±12 nm respectively) as 
checked by a traditional DLS apparatus (Brookhaven Instruments, ZetaPlus). The two smallest 
samples have been dispersed in water at a concentration of 0. 1 % w/w, whereas the largest ones 
have been prepared in water at a concentration of about 0.001% w/w. All the samples gave a 
sufficiently high scattering signal and for each sample, various sequences of 100 images were 
acquired at 7 different exposure times At ofabout0.5 ms, 1.5 ms, 5.5 ms, 18ms, 55 ms, 180 ms, 
and 550 ms. For each sequence, the ETDS was calculated. The optical background due to stray 
light contributions was removed from each image by subtracting the image average obtained 
over the set of images for each exposure time. In Fig. |2]we show the images obtained with the 
80 nm particles at four different exposure times. As it is apparent from the figures, as the expo- 
sure time increases the appearance of the images changes, and the short scale inhomogeneities 
are progressively smeared. 

In Fig. [3] the 5, (q,Af) spectrum for 80 nm particles is shown at different exposure times. 
This data show clearly that the spectra decrease as the exposure time is increased. For the 
longest exposure time the averaging of the fluctuations reduces the scattering signal towards 
the electronic background level. 

In Fig.|4]the 5, (q, At ) spectra are represented as a function of the exposure time at 4 different 
wave vectors. It's possible to observe the different decay times and the different amplitudes of 
the electronic background. Fitting curves obtained by using equations 4 and 5 are displayed, 
too. The fitting is carried out with three free parameters, namely, the time constant T (q), the 
product of the static scattered intensity times the instrument transfer function /(q) • T (q), and 
the electronic background of the system B(q). The agreement between data points and the 
corresponding fitting curves is very good. It's worth pointing out that, even if this fitting is made 
with such a limited number of data-points, each point is the output of a statistical analysis of 100 
images and several statistically independent samples, is the number of independent samples is 
about « w nq/q m in, where q m i„ is the minimum wave vector. This means that in the wave vector 




Fig. 2. Acquired images of nano-particles of 80 nm diameter at four different exposure 
times Af (a) 18 ms, b) 55 ms, c) 180 ms, d) 550 ms). The image size corresponds to 75 x 
75/lm 2 in real space. 



range of interest the number of statistical samples for each point in Fig. [4] varies from 3,000 to 
30,000. 

In Fig. [5] the time constants obtained by the proposed fitting procedure are plotted for the 
different samples as a function of the wave vector. Lines correspond to fitting of x values ac- 
cording to the formula %{Q) = 1/ (DQ 2 ). This fitting procedure has only D as an adjustable 
parameter; the particle diameter can be recovered from the Stokes-Einstein equation. Obtained 
diameters are 81 ±2 nm, 138±7 nm and 427±12 nm, in very good agreement with the sizes 
measured by standard DLS. 

5. Conclusions 

In the present work we have presented results from a dynamic near field scattering experiment, 
performed on colloidal particles by means of a simple set-up, consisting of a suitable modi- 
fied optical microscope and a laser beam as illuminating light. A new statistical algorithm is 
presented that allows us to measure the characteristic time constants of the colloidal system. 
The statistical procedure enables the camera to detect characteristic times much shorter than 
the camera delay times. Indeed, the only limitation is associated with the minimum exposure 




Fig. 3. Exposure Time Dependent Spectra S; (q, At) measured as a function of the scattering 
vector q, at different exposure times At for particles of 80 nm diameter. The central arrow 
indicates the increase of the exposure time, while the vertical lines indicate the lower and 
upper boundaries of the q range over which the fitting procedure has been performed. 
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Fig. 4. Exposure time dependent spectra 5/ (q, At) measured as a function of the exposure 
times At, at four different wave vectors q, for particles of 80 nm diameter. Lines correspond 
to fitting curves determined by Eq.[l0]and[TT]see the main text for details. 
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Fig. 5. Time constant X measured as a function of the wave vector q, for 3 different particle 
sizes. Straight lines in the log-log plot correspond to fitting curves with the power law 



T(Q) = l/(DQ 2 ). 



time of the camera and not by its delay time. This result opens the way to the investigation 
of ultra-fast dynamics such as molecular motions, capillary waves, nano-particles diffusion, 
and virus or biological macromolecules mobility. Until now these phenomena required sophis- 
ticated equipment, whereas the proposed method uses a simple modified microscope and a 
standard camera. 
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